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Abstract 

Inaccuracies  in  using  the  off-axis  tension  test  are 
caused  by  the  nonuniform  stresses  produced  by  the  end 
constraints.  The  purpose  of  this  -study-  i  s  to  show  that 
these  nonuniformities  can  be  virtually  eliminated  by 
adjusting  the  amount  of  tab  clamping  and  selectively 
locating  the  point  about  which  the  clamp  may  rotate.  A 
second  phase  of  the  study  is  to  conduct  a  limited 
verification  of  a  strain  energy  failure  criteria  that  can 
account  for  the  material  nonlinearity  found  in  composite)^/, 

P 

Determination  of  the  ideal  amount  of  clamping  and 
point  of  rotation  was  carried  out  using  a  linear  finite 
element  method.  The  nonlinear  behavior  of  the  specimen  was 
simulated  using  a  nonlinear  finite  element  method.  The 
failure  criteria  was  used  in  conjunction  with  the  nonlinear 
finite  element  algorithm  in  determining  failure. 
Experimentation  was  also  carried  out  to  verify  the  finite 
element  solutions. 

The  results  indicate  that  a  nearly  uniform  state  of 
stress  can  be  produced  with  ideal  tab  clamping  and 
rotation.  The  results  indicate  that  the  nonlinear  behavior 
of  the  specimen  was  accurately  predicted  by  the  nonlinear 
finite  element  method  but  that  the  failure  state  predicted 
by  the  failure  criteria  was  in  error.  Suggestions  are  made 


for  improving  the  accuracy  of  the  failure  criteria. 


IMPROVEMENT  OF  END  BOUNDARY  CONDITIONS 
FOR  OFF-AXIS  TENSION  SPECIMEN  USE 


CHAPTER  1 

INTRODUCTION 

1 .  I .  Background 

In  the  off-axis  tension  test  specimen,  unidirectional 
composite  material  fibers  are  oriented  at  an  angle  that  is 
neither  0°  nor  90°  to  the  specimen  axis.  If  the  load  is 
introduced  uniaxial ly  with  respect  to  the  specimen  axis, 
normal  and  shearing  stresses  will  develop  in  the  specimen 
relative  to  the  material  axes  1  and  2.  Thus,  it  would 
seem,  the  off-axis  tension  test  is  a  very  simple  means  of 
testing  fibei — reinforced  composites  under  a  biaxial  state 
of  stress.  Furthermore,  since  available  experimental  data 
confirm  the  assumption  that  unidirectional  composites  are 
orthotropic  i.e.  normal  and  shear  stresses  are  uncoupled, 
the  off-axis  test  would  appear  to  be  a  very  convenient 


method  for  testing  the  intralaminar  (in-plane)  shear 
propert i es  of  compos i tes . 

The  usefulness  of  the  off-axis  tension  test  Is 
somewhat  diminished  when  one  attempts  to  load  an  off-axis 
specimen  with  a  truly  uniform  axial  load.  Since  the  off- 
axis  specimen  is  anisotropic  with  respect  to  the  specimen 
axest  a  uniform  axial  stress*  o^,  causes  not  only  normal 
and  transverse  strains  but  also  shearing  strains  relative 
to  the  specimen  axes.  Consequently,  any  resistance  to 
deformation  due  to  end  constraints  such  as  tabs  or  clamps 
will  cause  the  axial  stress,  o  ,  to  become  nonuniform. 

This,  of  course,  also  causes  the  biaxial  stress  field  to 
become  nonuniform.  Thus,  any  test  results  obtained  using 
an  off-axis  specimen  will  most  likely  be  Influenced  by  this 
nonuniformity.  Though  the  influence  of  end  constraints  is 
present  to  some  extent  in  all  tensile  tests,  it  is 
magnified  in  the  case  of  off-axis  specimens  because  of  the 
rotational  tendency  caused  by  the  shear  coupling  effects 
(1).  Various  techniques  aimed  at  eliminating  this 
difficulty  have  been  suggested.  These  techniques  will  be 
discussed  in  the  next  few  paragraphs  in  approximately 
chronological  order. 

Tsai  (243)  used  off-axis  specimens  to  experimentally 
verify  his  strength  criteria.  The  criteria  is  based  on  the 
assumption  that  the  axial  stress,  be  uniform.  He  found 

that  specimens  of  uniform  cross  section  and  low  fiber 


orientation  angles  had  a  tendency  to  fail  under  the  rigid 
grips  used  to  apply  the  load.  To  Insure  that  failure 
occurred  in  the  test  section,  he  reduced  the  test  section 
in  a  "dog-bone”  fashion.  In  a  later  effort,  Lauraitis  (4) 
used  off-axis  specimens  to  verify  a  strength  criteria  based 
on  linear  elastic  fracture  mechanics.  For  the  same  reasons 
given  by  Tsai,  she  reduced  the  test  section  cross  section 
using  a  continuous  large  radius  of  curvature. 

Pagano  and  Hal  pin  (5)  presented  analytical  and 
experimental  evidence  of  the  adverse  effect  of  conventional 
(non-rotating)  clamping  devices  on  the  stress  field  in  an 
off-axis  specimen.  In  the  study,  they  presented  a  solution 
to  a  boundary  value  problem  that  approximated  the  off-axis 
specimen.  Experiments  with  nylon  reinforced  rubber  were 
also  carried  out  to  qualitatively  verify  the  analytical 
results.  The  study  showed  that  significant  nonuni  form  i n- 
piane  bending  stresses  are  produced  in  the  off-axis 
specimen  due  to  end  constraints.  After  allowing  limited 
rotat i on  of  the  c 1  amps ,  they  cone  I uded  that  the  gr i pp i ng 
restraint  was  the  dominant  factor  in  perturbing  the  stress 
field  and  that  clamp  rotation  was  of  little  importance. 

They  also  concluded  that  increasing  the  length  of  the 
specimen  would  produce  a  sufficiently  uniform  state  of 
stress  to  obtain  elastic  moduli.  To  obtain  strength  data, 
however,  they  pointed  out  that  some  modification  of  the 
ends  of  the  specimen  would  be  required. 


Rizzo  (1)  used  the  finite  element  method  to 
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investigate  the  effect  on  the  off-axis  specimen  stress 
field  of  completely  rigid  clamping  with  and  without  end 
rotation.  He  observed  a  marked  Improvement  in  stress 
uniformity  when  clamp  rotation  was  permitted.  He  also 
restated  the  observation  made  by  Hal  pin  and  Pagano  that  for 
long  specimens  ( 1 ength/width  >  10)  the  stress  field  at  the 
center  of  the  specimen  is  uniform  regardless  of  the  end 
clamping  arrangement.  Two  additional,  somewhat  subtle, 
observations  made  were  that  resistance  to  Poisson  type 
contraction  near  the  tabs  has  a  major  influence  on  the 
test's  inaccuracies  and  that  changing  the  pivot  point  may 
alter  the  above  situtation  significantly.  No  analytical  or 
experimental  data  was  given,  however,  to  support  these  last 
two  observations. 

Wu  and  Thomas  (6)  designed  a  test  fixture  which 
allowed  rotation  of  the  end  clamps  about  a  point  centered 
at  the  edge  of  the  clamp  area.  With  this  fixture,  they 
tested  15°  off-axis  specimens  with  length  to  width  ratios 
of  5,  4,  and  2.5.  They  acquired  strain  data  from  nine 
strain  gages  spaced  relatively  evenly  across  the  specimen. 
They  concluded  from  their  investigation  that  with  their 
rotating  clamp  fixture  and  within  the  length  to  width 
ratios  used,  the  state  of  strain  at  low  stress  levels 
appeared  to  be  relatively  uniform.  They  also  suggested 
that  putting  a  compliant  material  between  the  specimen  and 
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the  clamp  could  reduce  the  restraining  effect  of  the  rigid 
clamps  and  that  whether  or  not  the  fixture  could  provide 
the  desired  degree  of  uniformity  at  the  high  stress  levels 
required  in  strength  testing  remained  to  be  established. 

Richards  et.  al .  (7)  conducted  tests  using  45°  specimens 
with  three  inch  tapered  tabs,  length  to  width  ratios  of  12, 
and  loaded  without  any  end  rotation  allowed,  which  verified 
the  findings  of  Ha  I  pin  and  Pagano  concerning  the  effect  of 
large  length  to  width  ratios  on  uniformity. 

Cole  and  Pipes  (8)  examined  the  off-axis  tensile  test 
In  order  to  determine  its  utility  as  a  biaxial 
characterization  specimen  for  boron  epoxy.  Attempting  to 
insure  a  uniform  test  section  stress  field,  they  used  long 
specimens  with  the  fiber  orientation  of  highly  tapered  end 
tabs  the  same  as  the  specimen  fibers.  The  choice  of  tab 
fiber  orientation  was  supported  only  by  a  simple  Intuitive 
argument.  A  very  important  discovery  made  In  the  course  of 
their  study  was  the  fact  that  the  normal  and  shear 
responses  of  boron  epoxy  remain  uncoupled  through  the 
nonlinear  regime.  This  observation  prompted  the  advocation 
of  the  off-axis  specimen  for  use  in  determining  shear 
stress-strain  response.  Cham Is  and  Sinclair  (9)  analyzed 
the  10°  specimen  explicitly  for  that  purpose.  Their  test 
fixture  allowed  rotation  of  the  end-clamps  about  points 
past  the  ends  of  the  specimen.  Chang  et.  al.  (10)  tested  a 
modified  version  of  the  Chamis  and  Sinclair  fixture  which 


had  the  clamp  rotation  point  centered  at  the  tapered  edge 
of  the  tab.  This  modification  appeared  to  Improve  stress 
uniformity. 

Sandhu  and  SendeckyJ  (11)  conducted  an  extens 1 ve 
parametric  study  to  determine  a  tab  fiber  orientation  and 
tab  Inclination  angle  which  produced  the  most  uniform  state 
of  stress  In  off-axis  specimens  with  off-axis  angles 
ranging  from  4°  to  80°.  The  study  was  carried  out  with  a 
linear  finite  element  program  for  both  rotating  and  non- 
rotating  clamps.  They  found  that  matching  the  tab  fiber 
orientation  with  the  specimen  fiber  orientation,  as 
suggested  In  Ref.  (8),  did  In  fact  reduce  the  nonuniformity 
of  the  stress  field.  They  also  found  that,  using  an 
optimized  tab  fiber  angle  and  tab  Inclination  angle  along 
with  rotating  clamps,  a  nearly  uniform  state  of  stress 
could  be  attained.  It  is  important  to  note  that  the 
specimen  design  In  their  study  did  not  rely  on  large  length 
to  width  ratios  for  achieving  a  uniform  state  of  stress. 

To  verify  the  results  of  the  parametric  study,  experiments 
with  the  above  optimized  tab  parameters  were  carried  out. 
The  results  of  the  experiments  were  then  plotted  against 
the  response  predicted  by  a  nonlinear  finite  element 
program  incorporating  appropriate  boundary  conditions. 


1.2.  Purpose 


The  purpose  of  this  study  Is  to  show  that  a  state  of 
nearly  uniform  stress  can  be  produced  In  a  standard 
geometry  off-axis  specimen  by  adjusting  the  amount  of  tab 
clamping  and  selectively  locating  the  point  about  which  the 
clamp  will  rotate.  Since  In  the  analysis  phase  of  the 
study,  a  uniform  state  of  stress  In  the  entire  gage  length 
(area  between  the  tabs)  will  be  the  design  objective,  this 
study  will  result  In  a  test  method  that  will  be  useful  not 
on  1 y  for  obta 1 n i ng  el ast 1 c  modu 1 1  but  a  1  so  for  conduct  1 ng 
strength  tests. 

A  secondary  purpose  of  this  study  is  a  limited 
verification  of  the  strain  energy  failure  criteria 
presented  by  Sandhu  In  Ref.  <12  &  13). 

The  present  study  will  be  an  extension  of  the  work 
done  by  Sandhu  and  Sendecky j .  The  essential  difference 
between  the  two  studies  Is  the  parameters  that  are  allowed 
to  vary.  In  Sandhu  and  SendeckyJ's  work,  the  parameters 
were  fiber  orientation  angle  and  tab  Inclination  angle 
whereas  In  the  present  study  the  parameters  are  the  amount 
of  clamping  and  the  point  about  which  the  clamped  area  will 


rotate. 


CHAPTER  2 


I 


8 

i 


THEORY 


2.1.  I ntroduct  f  on 

In  this  section,  the  theory  used  to  describe  the 
response  of  composite  laminates,  as  It  relates  to  the  off- 
axis  specimen,  will  be  discussed.  Second,  a  formulation  of 
the  finite  element  methods  used  (linear  and  nonlinear)  will 
be  presented.  Finally,  the  failure  criteria  to  be  verified 
will  be  presented. 

2.2.  Macromechanical  Behavior  of  the  Off-Axis  Specimen 


c 


If  a  uniform  axial  stress,  o^,  is  applied  to  an  off- 
axis  specimen,  as  shown  in  Fig.  2.1,  normal  and  transverse 
stresses  will  develop  relative  to  the  material  axis  system. 
Recalling  from  elementary  mechanics  of  materials,  the 
transformation  equations  relating  the  stresses  in  an 
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X,Y  coordinate  system  to  the  stresses  In  a  1,2  coordinate 


system. 


or, 


J  \ 

m2 

n2 

2nm 
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-2nm 
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-nm 

nm 

m  -  n 

*xy 

S.  / 

{o}  = 

CT]{ 

o) 

(2.1) 


(2.2) 


where. 


n  =  sine  and  m  =  cose 
We  find  for  the  unlaxially  loaded  off-axis  specimen. 


(2.3) 


O  s  X  =0 
y  xy 


giving. 


1 


m^o. 


®2  =  "  ®x 


(2.4) 


T  j  2  =  nmo^ 


Obviously  then,  if  we  can  find  a  way  to  load  the  off- 
axis  specimen  with  a  truly  uniform  axial  stress,  we  will 
have  a  very  convenient  method  of  testing  unidirectional 
compos i te  mater  1  a  1 s  under  a  un I  form  biaxial  state  of 
stress . 

Note  however  (see  Fig.  2.1),  that  a  shearing  strain 
is  induced  relative  to  the  specimen  axes  when  is 
appi led.  This  Is  due  to  the  fact  that  the  off-axis 
specimen  is  anisotropic  with  respect  to  the  specimen  axes. 
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^66 
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It  should  be  pointed  out  that  the  stress-strain 
response  is  not  anisotropic  with  respect  to  the  material 
axes,  rather  it  is  orthotropic,  i.e.. 


Sji  Sj2 


S,2  ^22 


^66  ^12 


(2.6) 


Thus,  the  anisotropic  behavior  observed  relative  to 
the  specimen  axes  is  due  solely  to  the  transformation  of 
the  orthotropic  stress-strain  relations  through  an  angle 
9, 


[5]  =  [T]^rS][T] 


(2.7) 


where,  [S]  and  [T]  are  as  given  above.  All  the  terms  of 
the  [S]  and  [S]  matrices  are  given  explicitly  in 
appendix  A. 

In  order  to  physically  load  a  specimen  to  the  high 
load  levels  required  for  most  testing,  some  sort  of 
clamping  device  (such  as  thosse  found  in  an  Instron  or  MTS 
machine)  will  undoubtedly  be  used.  Furthermore,  the  ends 
of  the  specimen  may  be  reinforced  by  tabs  of  one  sort  or 
another.  If  the  loading  device  and  end-tabs  provide  any 
resistance  to  the  natural  deformations  associated  with  a 
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uniform  o^,  the  axial  and  biaxial  stress  fields  will  become 


nonun  I  form. 

The  shear  coupling  ratio  (14:22), 

(2.8) 

gives  a  measure  of  the  degree  to  which  end-rotation  will 
occur  as  a  function  of  off-axis  angle.  The  greater  the 
possible  rotation,  the  greater  the  effect  of  the  clamp.  If 


’’xy  ^16^^  11 


we  set  the  derivative  of  with  respect  to  the  off-axis 

X  X 


angle  equal  to  zero,  we  can  solve  for  the  angle  for  which 
shear  coupling  is  a  maximum. 

-  0  (2.9) 


dn^y/ae  =  <s,,s;6  -  Sissji )/(S, , )‘  = 


where  the  primes  denote  differentiation  with  respect  to  0. 
The  terms  in  Eq.  (2.9)  are. 


Su  =  +  m^n2(2Sj2+Sg^)  +  n'‘S22 


♦  ^2„2, 

1  ^ 

er  \  o _ 9_ 


d§jj/de  =  2mn^( 2522-25 J2-S66)  -  2m  n ( 25 j j -25 j 2-5^^ ) 


(2.10) 


^16  "  ^^12  ^66^  n^rn( 2522-25  j  2  5^^) 


dSi^/d0  -  m  (25jj  25 j 2  S^^)  6n^m  iS,j+S22  2bj2-Sgg) 


+  n  (2S22  25 j 2  Sgg) 


2  , 


51  nee  (Sj J )  is  positive  and  finite,  we  need  only  solve 


^11^16  ^16^1  1  °  (2.11) 


for  9.  This  equation  will  be  used  subsequently  to 
determine  the  off-axis  angle  that  maximizes  the  shear 
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V' 


coupling  ratio  for  the  material  used  in  this  study. 


The  work  done  by  Sandhu  (12)  and  Daniel  (15)  and 
reported  by  Whitney  et.  al .  (14:192)  indicate  that  the  off- 
axis  specimen  does  not  yield  as  much  of  the  shear  stress- 
strain  curve  as  does  the  (±45)  laminate  commonly  used  for 
obtaining  shear  stress-strain  curves  for  graphite/epoxy.  A 
possible  reason  for  this  phenomenon  is  that  the  method  of 
load  introduction  (rigid  clamps  with  little  or  no  rotation) 
cause  stress  concentrations  near  the  grips  that  in  turn 
cause  premature  failure.  In  this  study  it  will  be  observed 
if  the  near  elimination  of  such  stress  concentrations  will 
allow  the  off-axis  specimen  to  yield  a  shear  stress-strain 
curve  closer  to  that  of  the  (±45)  laminate. 

Refs.  (9)  and  (11)  point  out  that  the  fiber 
orientation  of  an  off-axis  specimen  used  for  obtaining 
shear  stress-strain  curves  must  be  chosen  such  that  the 
contribution  made  by  the  norma)  stresses  towards  failure  is 
a  minimum  or  equivalently,  the  contribution  made  by  shear 
stresses  is  a  maximum.  If  we  look  at  the  Tsai -Hi  11  failure 
cr i ter i a , 

(Oj/X)^  +  aj02/(X)^  =  1  (2.12) 

whe’"e  X,  Y,  and  5  are  the  normal,  transverse,  and  shear 
strengths  in  tension,  respectively,  we  see  that  the 
contribution  made  by  shear  is. 


■'.J 


2  2 

Tj2  =  -nmo^  ®2  ~  ^  (2.14) 

we  can  rewrite  Eq.  (2.12)  as, 

=  l/(m^/X^  +  (l/S^  -  l/X^)mV  +  n^/Y^)  (2.15) 

Thus,  substituting  Eq.  (2.15)  into  Eq.  (2.13)  we  obtain, 

=  1/((S  cot8/X)^  +  (1  -  (S/X)^)  +  (S  tane/Y)^)  (2.16) 

Now  Zg  attains  a  maximum  value  when, 

dZ^/da  =0  (2.17) 

which,  on  simplification  leads  to, 

tane'  =  \/(Y/X)  (2.18) 

where  a'  is  the  off-axis  angle  for  which  Z^  attains  a 
maximum.  This  equation  will  be  used  subsequently  to 
determine  the  value  of  9'  for  the  material  used  in  this 
study. 

It  is  important  to  note  that  all  of  the  derivations  to 
this  point  have  been  based  on  the  assumption  that  the 
material  under  consideration  remains  linearly  elastic  to 
failure.  This  is  not  the  case,  however.  In  reality,  the 
shear  response  of  graph i te/epoxy ,  the  material  system  used 
in  this  study,  is  highly  nonlinear.  Ref.  (11)  showed, 
using  a  fully  nonlinear  finite  element  technique,  that  the 
value  of  9'  corresponding  to  maximum  shear  contribution  to 
failure  was  13°  whereas,  using  the  linear  analysis  shown 
above,  the  value  of  9'  was  11.38°.  This  adjustment  must  be 


taken  into  consideration  when  using  Eq.  (2.18). 
2.3.  Linear  Finite  Element  Theory 


The  1  f ne^ar  finite  element  program  used  in  this  study 
is  based  on  the  constant  strain  triangle  (see  Fig.  2.2  for 
coordinates).  The  assumed  element  displacement  field  is 


1  i near  i n  x  and  y : 


The  strain  are  expressed  by. 


where  {d}  is  the  nodal  displacement  vector  and. 


(2.19) 


(2.20) 


[B]  =  1/2A  0  (x^-Xp)  0  (X2~^l^ 

(X3-X2)  (y2-y3>  ^y3~y0  (yi“y2> 


(2.21  ) 


where  A  is  the  area  of  the  triangle  and  x,  and  y.  are  the 


X 


9^. 


and  y  coord } nates  of  the  5th  node. 

For  an  orthotropic  lamina  whose  principal  direction  is 
oriented  at  an  angle  Q  to  the  X  and  Y  directions,  the 
stresses  are  related  to  the  strains  by. 


Now  allowing  for  n  layers  of  orthotropic  material 
through  the  thickness  we  can  express  the  equivalent  element 
stiffness  matrix  as, 

n 

S  [B]'^[Q],[B]t,  (2.Z6) 

I  =  I 

where  i  denotes  the  ith  layer.  Thus,  loads  and 
displacements  are  related,  on  an  element  basis  by. 
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(2.27) 


I 


I 

i 


? 


in  =  C»<eq3td} 

Note  that  even  though  multiple  layers  have  been 
allowed  through  the  thickness  of  the  element,  the  layers 
have  common  nodes  and  thus  no  strain  or  displacement 
variation  Is  being  permitted  through  the  thickness. 


I 

I  ’■ 


i 


2.4.  Nonlinear  Finite  Element  Theory 

The  effect  of  material  nonlinearity  on  the  uniformity 
of  the  off-axis  specimen  stress  field  Is  analyzed  using  a 
nonlinear  finite  element  technique  based  on  the  constant 
strain  element  described  in  the  preceding  section.  A 
complete  development  of  the  technique  Is  given  in  Ref. 
(12).  In  this  technique. 


•V 


u 


d{f}  »  Ck(c)gq]d{d) 


(2.28) 


where  d{f},  d{d} ,  and  Ck(e)^q]  are  the  increments  of  load 


and  displacement  and  the  current  stiffness  matrix, 
respectively.  As  shown,  the  stiffness  matrix  depends  on 
the  current  strain  level.  This  dependence  comes  about 


through  the  [Q]  matrix.  Recall  the  equation  for  [k^^]. 


(Eq.  2.26), 


[Xeql  ■  *  S 

i  =  l 


(2.29) 


One  can  see  that  [^<gq]  is  calculated  using  [Q]  ^  (See 


appendix  A  for  the  terms  of  [0]).  If  the  material  being 
analyzed  exhibits  material  nonlinearity,  the  lamina 
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material  properties  (Ej,  £2*  Gi2»  '*12^  will  vary  with 
strain.  It  follows  then,  since  [Q].  is  calculated  from 
these  properties,  that  [Q] j  and  hence  [kgq]  will  vary  with 
strain  as  well. 


i 


•  C 


fr 

• 


■r, 


►5 


As  Implied,  the  material  property  curves  for  the 
material (s)  being  analyzed  are  obtained  from  the 
experimental  stress-strain  curves  of  unidirectional  lamina. 
The  curves,  and  corresponding  experimental  tests,  required 


5"* 

•k:- 

for  each  material 

1  ,  are  shown 

1 n  Table  2.1. 

Curve 

Experiment 

vs . 

0°  tension 

vs . 

0®  compression 

. 

vs . 

®2 

90°  tension 

b 

‘’z 

vs . 

£2 

90°  compression 

^12 

vs . 

^12 

t45°  tension 

1 

j* 

'’12 

vs . 

0°  tension 

''12 

vs . 

e . 

i 

0°  compression 

Table  2.1  Stress-Strain  Curves  and 
Corresponding  Experimental  Tests 


To  make  the  exper  1  ment.3 1  curves  usable  by  the  computer 
program,  they  are  represented  analytically  by  cubic  spline 
Interpolation  functions.  The  material  properties  curves 
are  thus  represented  as  functions  of  strain  by  simple 
polynomials.  With  the  stress-strain  curves  In  this  form. 
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the  tangent  moduli  £2*  ^12^  can  be  read'ly  obtainedt 

as  functions  of  strain,  by  differentiating  the  appropriate 
cubic  spline  functions  with  respect  to  strain. 

Obviously,  the  load  cannot  be  applied  in  a  truly 
continuous  manner  using  a  computer  program.  As  indicated 
in  Eq.  (2.28),  It  Is  applied  Incrementally.  This  presents 
a  slight  dilemma,  however,  since  may  vary  during  a 

load  increment.  To  overcome  this  difficulty,  a  predictoi — 
corrector  Iterative  procedure  is  utilized.  In  this 
procedure  is  initially  calculated  using  material 

properties  that  correspond  to  the  state  of  strain  existing 
at  the  end  of  the  previous  load  increment.  For  the  very 
first  increment,  engineering  1  inear  elastic  constants  are 
used. 

Having  calculated  Ck^^],  corresponding  to  the  end  of 
the  previous  load  increment,  a  new  increment  of  the  load  is 
applied  and  an  increment  of  displacement,  d{d} ,  is 
calculated  using  Eq.  (2.28).  Note  that  when  a  given  load 
increment  is  first  applied,  Eq.  (2.28)  will  take  the  form, 

d{f}n^l  =  (2.30) 

where  n  denotes  the  nth  load  increment. 

d{d}  can  then  be  used  to  calculate  an  increment  of 
strain  using  a  slightly  modified  version  of  Eq.  (2.20), 


d{e)  =  [B]d{d) 


(2.31  ) 


The  increment  of  stress  can  be  calculated  from, 

d{o),  =  [Q],d{c}  (2.32) 

were  I  denotes  the  1th  layer  as  defined  In  section  2.3. 

Notice  that  In  this  last  equation  It  has  been  assumed 

that 

a.  The  Increment  of  stress  depends  upon  the  strain 
state  and  the  Increment  of  strain;  and 

b.  The  Increment  of  strain  Is  proportional  to  the 
Increment  of  stress. 

Note  that  the  stress  and  strain  Increments  In  Eqs. 
(2.31)  and  (2.32)  are  relative  to  the  specimen  axes  X  and 
Y.  To  obtain  the  stress  and  strain  increments  relative  to 
the  material  axes,  1  and  2,  Eq.  (2.1)  must  be  used, 

d{a)  =  [T]d{a)  (2.33) 

d{e}  =  CT]d{e)  (2.34) 

Having  computed  the  Increments  of  stress  and  strain, 
new  levels  of  stress  and  strain  are  determined  by  adding 
the  new  Increments  to  the  levels  existing  at  the  end  of  the 
previous  increment. 

A  mean  level  of  strain  Is  then  computed  by  averaging 
the  new  level  of  strain  with  the  level  of  strain  existing 
at  the  end  of  the  previous  load  increment.  These  "average" 
strains  are  then  used  to  determine  a  new  set  of  material 
properties  since,  as  mentioned  above,  they  are  readily 
obtainable  as  functions  of  strain  through  the  cubic  spline 


functions.  This  concept  will  be  discussed  graphically  in  a 


subsequent  paragraph. 

With  the  new  elastic  properties  in  hand,  [Q]  and  hence 
Ckgq]  are  recalculated.  The  toad  increment  Is  then 
reapplied  and  the  Increment  of  displacement  and  strain,  is 
recalculated.  Note  that  the  same  load  increment  has  been 
reapplied,  not  an  additional  one.  This  procedure  Is 
continued,  for  a  given  load  Increment,  until  the  ratio  of 
the  change  of  the  strain  Increment  to  the  strain  increment 
is  less  that  0.001,  l.e., 

~  ^  0-001  (2.35) 

where  n  denotes  the  nth  reappl Icatlon  of  a  given  load 
increment.  When  Eq.  (2.35)  Is  satisfied,  a  new  load 
increment  is  applied.  The  repetitive  use  of  the  procedure 
outlined  above  generates  the  predicted  stress-strain 
response  of  the  laminate  under  consideration. 

In  Fig.  2.3  an  attempt  has  been  made  to  graphically 
display  how  the  above  procedure  is  carried  out.  For 
clarity  a  simple  case  of  one-dimensional  loading  in  shear 
is  shown.  As  shown,  an  initial  strain  Yq  's  assumed  to 
exist  due  to  a  previous  load  Increment.  An  increment  of 
the  load  is  then  applied  and  the  resulting  strain,  Yj,  Is 
calculated  based  on  the  tangent  modulus  at  point  0.  A  new 
modulus  is  now  determined,  corresponding  to  the  average  of 
Yj  and  Yg  (point  A).  The  load  Is  then  reapplied  and  the 


is  calculated  based  on  the  modulus  at  A. 


strain, 

Again,  the  modulus  corresponding  to  the  average  of  Y2  and 
Yq  (point  B)  is  determined  and  the  load  reapplied,  giving 
rise  to  Yj-  Por  this  simple  case  the  procedure  is  repeated 
for  the  given  load  increment  until, 

^^n+l  ■  ^  (2.36) 

For  graphite/epoxy  and  glass/epoxy,  the  materials  used 
in  this  study,  nonlinearity  is  most  pronounced  in  the  shear 
direction  (relative  to  the  material  axes)  while  in  the 
normal  and  transverse  directions  the  stress-strain  response 
is  nearly  linear  (See  appendix  B  for  experimental  stress- 
strain  curves).  While  it  might  be  reasonable  to  assume  the 
normal  and  transverse  responses  to  be  linear,  such  an 
assumption  is  not  made  in  this  study.  Any  nonlinearity 
that  exists  in  any  of  the  material  property  curves  will  be 
accounted  for.  An  assumption  that  is  made,  however,  is 
that  nowhere  in  the  off-axis  specimen  will  the  strains  be 
of  a  large  enough  negative  value  to  justify  using  true 
compression  curves  as  indicated  in  Table  2.1.  Instead,  the 
curves  obtained  for  tension  will  be  used.  This  assumption 
was  found  to  be  valid  when  the  specimen  was  analyzed. 

In  the  procedure  Just  described,  the  lamina  biaxial 
strains  (ej  and  £2)  are  modified  before  being  used  to 
determine  the  elastic  constants  from  the  cubic  spline 
stress-strain  curves.  This  modification  is  required  to 
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allow  for  the  simultaneous  existence  of  longitudinal  and 
transverse  stresses  In  the  lamina  whereas  In  the 


experimental  data  only  one  component  of  stress  Is  present. 
Consider  the  differential  element  of  unidirectional 
material  shown  In  Fig.  2.4  (from  Ref.  (12)),  subject  to 
biaxial  stresses  and  02*  Recalling  that  we  have  assumed 
an  Increment  of  strain  Is  proportional  to  an  increment  of 
stress  we  can  write, 

d£  j  =  do  j  /E  j  ~  v  j  j  ~  ^  2  ^  (2.37) 

d62  ~  d02/E 2  ~  '*21  1  ^^2  ~  ^®2^^2 ^  ^ ~'*2 1  ^ 1  ^ ^  ^  (2.38) 

Eqs.  (2.37)  and  (2.38)  reveal  that  It  would  be 
erroneous  to  use  dtj  or  de2*  corresponding  to  the  biaxial 
stress  state  do^,  do2,  associated  with  an  off-axis 
specimen,  to  determine  Ej  or  E2  from  stress  strain  curves 
obtained  under  simple  loading  conditions.  For  example,  de2 
of  equation  (2.38)  corresponds  to  the  curve  ON  (Fig.  2.4) 
on  the  plane  OEHG  while  the  simple  stress  strain  curve  OM 
lies  on  the  plane  OEDC.  Since  stress-strain  data  similar  to 
ON  is  not  available,  we  have  assumed  that  simple  equivalent 
strain  Increments  can  be  computed  from  the  following 
express i ons : 


de  1 

eq 

=  d<j,/Ej 

=  de 1 / ( 1 -w , 2 (da2/do j ) ) 

(2.39) 

de2 

eq 

=  da2/E2 

=  de2/( l-v^j (daj/do2) ) 

(2.40) 
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The  error  that  would  result  from  using  de^  instead  of 


^^ijeq’  written  from  Eqs.  (2-39)  and  (2.40)  as, 

^‘^^ijeq  "  eq  =  (2-41) 

^‘^^2|eq  -  ^"2>/'^=2|eq  =  ^21  '2.42) 

From  Eq.  (2.4)  we  can  write,  for  the  off-axis  specimen, 

dOj  =  m^do^  (2.43) 

da  2  =  (2.44) 

which  upon  substituting  into  Eqs.  (2.41)  and  (2,42)  gives, 

Vj2  (clo2^^®l  ^  (2.45) 

V2j(doj/do2)  -  (2.46) 

To  get  a  feel  for  the  errors  that  would  result  in 


neglecting  this  correction  let  us  look  at  the  errors  that 
would  arise  in  a  14®  off-axis  specimen.  For  the  material 
properties  used  in  this  study  we  get, 

Vj2('tan29)  =  0.0174  =  1.74  %  error 
v^jlcot^s)  =  0.3514  =  35.14  %  error 

2.5.  Failure  Criteria 

The  incremental  loading  technique  described  in  the 
preceding  section  is  a  finite  process,  culminating  in  the 
failure  of  the  lamina.  To  determine  the  conditions  which 


result  in  failure,  various  criteria  have  been  proposed. 
These  criteria  assume  linear  material  behavior  to  failure 


As  mentioned  earlier,  however,  this  is  not  an  accurate 
assumption  for  graphIte/epoxy.  Consequently  a  failure 
criteria  has  been  developed  by  Sandhu  (11  &  12)  which  can 
account  for  nonlinear  material  behavior. 

Following  the  development  found  in  Ref.  (12),  a  scalar 
function,  f,  defining  the  failure  condition  of  materials 
exhibiting  nonlinear  behavior  can  be  written  as, 

f(6,o,K)  =  1  (2.47) 

where  e  and  o  are  the  stress  and  strain  states  and  K 
represents  the  material  character i st ics . 

The  explicit  form  of  Eq.  (2.43)  proposed  in  Ref.  (12) 
uses  the  scalar  strain  energy  to  determine  the  effect  of 
both  stress  and  strain  states  on  the  material  behavior. 

The  fact  that  this  criteria  is  a  function  of  both  stress 
and  strain  states,  helps  to  account  for  the  nonlinear  shear 
strains  that  can  occur  in  graph! te/ epoxy  prior  to  failure. 

Assuming  strain  energies  are  independent  parameters, 
the  failure  criterion  for  orthotropic  materials  may  be 
expressed  for  the  plane  stress  condition  as. 


where , 


(Kj/Kj)'"  +  (K2/K2)'"  +  (Kg,/Kg,)'"  =  I 


(i  =  1,2,6) 


(2.48) 


(2.49) 


s 
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and,  using  the  results  of  tests  under  simple  load 


condi t i ons , 


(i  =  1,2,6) 


(2.50) 


where  e.  are  the  current  strain  components,  are  the 

ultimate  normal  and  shear  strains  from  the  simple  loading 
tests  and  m  is  a  parameter  defining  the  shape  of  the 
failure  surface  in  the  strain  energy  space.  The  m 
parameter  i s  to  be  chosen  so  that  exper imenta 1 -ana  1 yt i ca 1 
correlation  is  best. 

An  analytical  comparison  of  various  failure  criteria. 
Including  the  present  one,  is  shown  in  Fig.  2.5  for 
boron/epoxy  (from  Ref.  (12)).  The  criteria  is  shown  for 
three  values  of  m,  namely  m  =  1/2,  1,  and  2.  One  can  see 
that,  in  the  Oj,02  the  failure  envelope  is  a  smooth 

continuous  curve  for  m  =  1.  For  m  >  1  the  curve  approaches 
a  discontinuity  as  is  found  in  the  maximum  stress  failure 
criteria.  It  may  be  true  that  some  value  of  m  <  1  would 
result  in  a  more  accurate  failure  criteria  than  m  =  1,  but 
since  there  is  no  biaxial  strain  energy  failure  data 
available  to  fix  m  to  a  specific  value  it  will  be  taken  to 
be  unity  as  was  done  by  Sandhu  in  Ref.  (12).  This  reduces 
the  criteria  to  a  simple  linear  relationship  of  strain 


energy  ratios. 
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CHAPTER  3 


ANALYSIS 


3.1.  I ntroduct 1  on 


The  analysis  of  the  off-axis  specimen  was  carried  out 


in  two  parts.  A  parametric  study  was  carried  out  to 


determine  the  amount  of  tab  clamping  and  point  of  rotation 


of  the  clamped  area  that  yielded  the  most  uniform  state  of 


stress  in  the  specimen  test  section.  Then,  a  nonlinear 


finite  element  analysis  was  carried  out  to  assess  the 


effect  of  material  nonlinearity  on  stress  field  uniformity. 


First,  however,  a  discussion  of  the  specimen  geometry  and 


finite  element  modeling  is  in  order. 


3.2.  Specimen  Geometry 


The  specimen  geometry  under  consideration  is  shown  in 


Fig.  3.1.  The  coupon  is  fabricated  from  16-ply 


3.1  Specimen  Geometry 


graphite/epoxy  and  is  tabbed  with  bi-directional 
glass/epoxy  tabs  (See  appendix  C  for  experimental  material 
properties).  Note  that  the  tab  fibers  are  aligned  parallel 
and  perpendicular  to  the  coupon  fibers  as  suggested  In 
Refs .  8  and  1 1 . 

Only  two  off-axis  angles,  namely  10°  and  14°,  were 
analyzed.  The  study  was  limited  to  only  two  angles  to 
conserve  computer  resources  and  still  allow  for  a 
comparison  of  the  optimized  boundary  conditions  at 
different  angles. 

The  first  reason  for  choosing  these  particular  angles 
is  that  they  bound  the  peak  values  of  the  shear  coupling 


ratio,  Hv-w*  From  Eq.  (2.11),  shown  here  again  for  clarity, 
xy 


^11^16'  ^16^11“  ° 


(2.11) 


we  find,  upon  solving  for  9,  that  for  the  experimental 


material  properties  under  consideration,  attains  its 

xy 


maximum  value  at  an  off-axis  angle  of  12-47  (Note  that  Eq. 
(2.11)  was  solved  numerically  using  a  simple  secant 
iterative  algorithm).  Furthermore,  looking  at  Fig.  3.2  one 
can  see  that  a  definite  peak  occurs  in  the  shear  coupling 
ratio  at  this  angle.  Thus,  by  optimizing  the  specimen  at 
these  two  angles,  close  bounding  to  a  worst  case  condition 
for  shear  coupling  was  obtained. 

A  second  reason  for  choosing  these  angles  Is  that  the 
shear  contribution  to  failure  is  maximized  in  this  range  of 
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)hear  Coupling  Ratio 


off-axis  angles.  Equation  (2.18) 


n'. 
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S 


i 
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tan  e'  =y  Y/X  (2. 18) 

yields  a  value  for  8'  of  8.5®,  using  experimentally 
determined  values  of  Y  and  X  in  tension.  Adjusting  this 
va ) ue  for  mater i a  1  non  1 i near  1 ty  by  the  same  amount 
indicated  by  Ref.  (11)  (See  the  end  of  section  2.2),  9' 
becomes  10.1®,  which  falls  into  the  10®  to  14®  range. 

Note  in  Fig.  3.1  that  the  specimen  length,  L,  is 
different  for  each  off-axis  angle.  The  length  is 
determined  from  the  requirement  that  a  fiber  passing 
through  the  center  of  the  specimen,  begin  and  end 
approx i mate  I y  one  I nch  from  each  tab .  This  feature  i s 
important  since  past  experience  indicates  that  fracture 
will  occur  along  a  fiber /matrix  interface.  The  length 
requirement  insures  that  such  a  fracture  can  occur  without 
constraint  from  the  tabs. 

A  specimen  width  of  one  inch  was  chosen  for  two 
reasons.  The  first  is  that  one  inch  is  a  standard  width 
for  composite  tensile  tests.  The  second  reason  is  that  it 
reduces  the  length  to  width  ratio  to  6  and  8  for  10°  and 
14°,  respectively.  By  designing  a  specimen  with  such  a 
small  length  to  width  ratio,  it  is  possible  to  observe 
whether  or  not  the  need  for  large  length  to  width  ratios 
can  be  eliminated  by  proper  application  of  the  load. 

The  remaining  features  of  the  specimen  geometry  are 
essentially  standard  and  self-explanatory.  One  final  note, 
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however,  fs  that  the  tab  taper  angle  was  adjusted  to  14° 
for  no  other  reason  than  to  allow  for  convenient  finite 
element  dimensioning.  Typical  values  for  the  taper  angle 
are  15°  and  20°. 

3.3.  Finite  El ement  Mode  1 1 ng 

A  good  finite  element  model  of  any  structure  should 
accurately  predict  the  displacement  field  of  the  actual 
physical  structure  under  consideration,  when  a  load  is 
applied  to  It.  To  obtain  such  a  model,  the  analyst  must 
accurately  model  the  size,  the  shape,  and  the  loading  (or 
displacement)  boundary  conditions  of  the  actual  structure. 
Furthermore,  he  must  provide  for  adequate  refinement  of  the 
element  mesh  In  areas  where  gradients  In  the  displacement 
field  are  known  or  expected  to  be  high.  In  describing  the 
finite  element  modeling  carried  out  In  the  present  study, 
each  of  the  characteristics  Just  mentioned  will  be 
dl scussed. 

Sizing  and  shaping  of  the  finite  element  model  In  the 
width  direction  was  obviously  trivial,  given  the  simple 
rectangular  pi  an form  of  the  off-axis  specimen.  Sizing  in 
the  thickness  direction  was  carried  out  by  specifying  the 
average  thickness  of  each  layer  of  material  In  a  given 
element,  1 .e. ,  tj  of  Eq.  (2.26).  In  the  portion  of  the 
specimen  where  there  Is  only  graph  1 1 e/epoxy  through  the 


thickness,  an  element  had  one  layer  with  a  thickness  of 
0.083  in.  In  the  tabbed  portion  of  the  specimen  where 
there  is  both  graphite/epoxy  and  glass/epoxy  through  the 
thickness,  but  no  clamp,  an  element  had  one  layer  of 
graphite/epoxy,  with  a  thickness  of  0.083  in.,  and  one 
layer  of  glass/epoxy  with  a  thickness  of  0.094  in.  In  the 
tapered  portion  of  the  tab,  the  thickness  of  the 
glass/epoxy  layer  was  reduced  in  a  stepwise  manner  from 
0.094  in.  to  0.000  in.  In  the  clamped  portion  of  the  tab, 
a  third  layer  of  steel,  1.5  in.  thick,  was  added  to  each 
element.  One  may  have  noticed  that  the  thickness  of  the 
glass/epoxy  has  just  been  given  as  0.094  in.  when  in  Fig. 
3.1  it  is  given  as  0.124  in.  This  is  due  to  a  slight 
modification  of  the  specimen  that  was  required  during 
experimentation.  This  modification  will  be  discussed  in 
chapter  4.  In  an  effort  to  keep  the  experimental  and 
analytical  specimens  as  si mi  liar  as  possible,  this 
modification  was  accounted  for  in  the  finite  element  model. 

The  experimental  test  fixture  used  in  this  study  to 
test  off-axis  specimens,  applies  the  load  with  rigid  Jaws 
that  are  free  to  rotate  about  a  pin  and  bearing  embedded  in 
the  Jaw.  This,  in  effect,  forces  all  points  in  the  clamped 
area  of  the  tab  to  rotate  concentrically  about  the  pin.  By 
adjusting  the  location  of  the  specimen  relative  to  the 
Jaws,  the  position  of  the  pin,  relative  to  the  clamped 
area,  can  be  adjusted.  A  detailed  description  of  the  test 


fixture  will  be  presented  in  chapter  4. 


To  model  this  type  of  rigid  rotating  boundary 
condition,  a  third  layer  of  steel  is  added  in  the  clamped 
area,  as  mentioned  above.  The  effect,  of  this  essentially 
rigid  layer,  is  to  constrain  the  clamped  portion  of  the  tab 
to  remain  rectangular  and  rotate  about  the  pin,  as  is  true 
for  the  actual  test  fixture.  Thus,  to  apply  a  load,  only 
the  displacements  at  the  two  nodes  corresponding  to  the  pin 
locations  (one  at  each  end  of  the  specimen)  need  to  be 
specified.  At  one  end  of  the  specimen,  the  displacements 
at  the  node  corresponding  to  the  pin  location,  are  fixed  at 
zero  (u  =  V  =  0)  At  the  opposite  end  of  the  specimen  a 
finite  displacement  is  applied,  at  the  node  corresponding 
to  the  other  pin  location,  in  the  X  direction  while  the  Y 
displacement  is  fixed  at  zero  (u^O,  vsO).  fig.  3.3 
shows  graphically  how  the  boundary  conditions  are  applied. 

The  work  of  Ref  (16)  clearly  shows  that  large 
gradients  can  exist  in  the  stress  field  of  the  off-axis 
specimen,  immediately  adjacent  to  the  tabs.  Though  the 
purpose  of  this  study  is  to  apply  the  load  in  such  a  mahher 
as  to  eliminate  these  gradients,  provisions  are  made  to 
accommodate  a  worst  case  condition  where  such  gradients  do 
exist.  As  can  be  seen  in  Fig.  3.3,  the  mesh  is  highly 
refined  near  the  edge  of  the  tabs.  The  mesh  is  also 
refined  in  the  tapered  portion  of  the  tab  so  that  a  smooth 
transition  from  glass/epoxy  and  graphite/epoxy  to  Just 
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graph f te/epoxy  Is  obtained. 

To  Insure  that  the  final  finite  element  model  design 
would  yield  results  that  had  converged  to  a  solution,  a 
convergence  study  was  carried  out  with  the  14°  specimen  In 
arriving  at  the  final  model.  The  three  models  used  are 
shown  In  Fig.  3.4.  As  a  measure  of  convergence,  axial 
stress  along  the  centerline  of  the  models  Is  plotted  In 
Fig.  3.5.  One  can  see  that  the  final  model  has  converged 
to  a  solution  since  at  all  points  along  the  curve,  the 
difference  between  the  medium  and  fine  models  Is  less  that 
5.0  %.  It  should  be  noted  that  the  quadrilateral  elements 
shown  In  all  of  the  finite  element  models  are  subdivided 
Into  four  constant  strain  triangles.  The  triangles  are 
formed  by  adding  a  node  at  the  center  of  each 
quadrilateral.  The  Interior  node  Is  subsequently  removed  by 
stat 1 c  condensat I  on . 

3.4.  Parametric  Analysis 

A  parametric  study  was  carried  out  using  the  finite 
element  technique  of  section  2.3  to  determine  the  amount  of 
tab  clamping  and  point  of  rotation  of  the  clamped  area 
which  produce  the  most  uniform  state  of  stress  In  the  gage 
length  (area  between  the  tabs)  of  the  off-axis  specimen. 

The  clamped  area  was  assumed  to  extend  across  the 

entire  width  of  the  tab  for  a  distance  D_,  measured  from 
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Fine  Mesh  1840  nodes 

Fig.  3-4  Finite  Element  Models  used  for  Convergence  3520  elements 


Irlot  of  Axial  Stress  Along  Specimen 
Centerline  to  Show  Convergence 


the  end  of  the  specimen,  as  shown  In  Fig.  3.6. 

The  pin  (or  point  of  rotation  of  the  clamped  area)  was 
assumed  to  be  located  along  the  center! ine  of  the  clamped 
area,  a  distance  Dp  from  the  end  of  the  specimen.  Note 
that  the  pin  can  only  be  located  at  the  boundaries  of,  or 
within,  the  clamped  area.  This  is  due  to  the  physical 
1  imitations  of  the  experimental  test  fixture. 

For  the  parametric  analysis,  the  clamping  distance, 

D^,  was  allowed  to  vary  from  I  in.  to  1  3/4  in.  For  a 
given  value  of  0^,  Dp  was  allowed  to  vary  from  0  to  D^,. 


1”  <  0^  ^ 

0”  <  D  <  D 
P  c 

where  D^  and  Dp  varied  in  1/8  in.  increments  since  in  the 
finite  element  model  the  clamp  must  terminate  at  an  element 
edge  and  the  pin  must  be  located  at  a  node. 

D^  was  limited  to  a  minimum  of  1  in.  due  to  the 
physical  limitations  of  the  fixture  and  specimen.  It  was 
felt  that  clamping  less  than  1  in.,  would  very  probably 
result  in  the  specimen  pulling  out  of  the  clamps  during 
test  I ng. 

As  a  measure  of  stress  field  uniformity,  the 
difference  between  the  maximum  and  minimum  axial  stress 
occur ing  anywhere  in  the  gage  length,  Ao^,  was  plotted 
versus  pin  position  for  a  given  amount  of  clamping.  The 
ideal  amount  of  clamping  and  position  of  the  pin  correspond 


vs 


to  the  values  of  D  and  0_  for  which  Ao  is  a  minimum.  The 

CP  ^ 

curves  for  the  10°  specimen  are  presented  in  Fig.  3.7. 

Those  for  the  14°  specimen  are  presented  in  Fig.  3.8.  Note 
that  the  results  for  each  pin  position  were  normalized  to  a 
center  section  axial  stress  of  30,000  psi  to  show  relative 
differences  in  Ao^. 

One  might  ask  whether  forcing  Aa^  to  go  to  zero  will 
also  cause  the  other  components  of  stress  (o  and  t  )  to 

y  Xy 

also  go  to  zero.  If  Ao^  is  zero  in  the  gage  length  then  o 
will  be  a  constant  in  the  gage  length.  From  the  x 
direction  plane  stress  equation  of  equilibrium. 


30j^/3x  +  atj^y/ay  =  0 


(3.1 ) 


we  can  see  that  if  is  a  constant,  can  only  change  in 

the  X  direction.  This  forces  t  to  be  zero  due  to  the 

xy 

stress  free  edges.  From  the  y  direction  equation  of 
equi 1 i br i um. 


30y/3y  +  3i^y/ax 


(3.2) 


one  can  see  that  if  t  is  zero,  o  must  also  be  zero  due 

xy  y 

to  the  stress  free  edges.  Thus,  by  driving  Ao  to  zero,  a 

X  y 

and  T  are  also  driven  to  zero, 
xy 

One  can  see  that  for  each  amount  of  clamping  there  is 


a  minimum  value  of  Figure  3.9  shows  this  minimum 

value  of  Ao^  plotted  versus  amount  of  clamping  (D^,),  for 
the  10°  and  14°  models,  respectively. 


Differential  Stress,  Act  (psi) 
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3.5.  Non-L  f  near  Ana  I ys  f  s 


. 
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To  assess  the  effect  of  mater  fa  1  nonf  fnearity  on  the 
off-axis  specimen  stress  field,  the  specimen  was  analyzed  a 
second  time  using  the  nonlinear  finite  element  technique 
presented  In  section  2.4.  For  this  analysis,  only  the 
1 dea 1  boundary  cond 1 t 1 on  for  each  of f-ax 1 s  ang 1 e , 
determined  In  the  preceding  section,  were  considered. 

From  the  linear  parametric  study  Just  discussed,  one 
can  see  that  the  ideal  amount  of  clamping  for  both  of f-ax Is 
angles  was  1  In.  The  ideal  pin  positions,  however,  were 
not  the  same.  For  the  10®  case,  the  Ideal  pin  location  was 
at  Op  ®  1/4  whereas  for  the  14®  case  the  Ideal  pin  position 

was  at  D  =0. 

P 

Note  that  the  Ideal  boundary  conditions  were 

determined  using  a  parametric  analysis  that  assumed  linear 

material  properties.  The  author  recognizes  that  due  to  the 

nonlinear  behavior  of  the  graph  1 te/ epoxy  off-axis  specimen, 

this  assumption  may  result  In  a  certain  amount  of  error  in 

determining  Ideal  values  of  D  and  D„.  The  work  of  Ref. 

c  p 

(11)  indicates  that  the  uniformity  of  the  stress  field  is 
degraded  by  material  nonlinearity.  That  study  also 
indicates,  however,  that  the  optimized  boundary  conditions 
arrived  at  assuming  material  linearity  are  a  good 
approximation  to  the  more  accurate  boundary  conditions  that 
are  obtained  if  the  nonlinear  technique  Is  used  for  the 


parametric  analysis.  Considering  that  the  ratio  of 


computer  cpu  seconds  for  a  nonlinear  run  to  cpu  seconds  for 
a  linear  run  is  approximately  150  (3000  sec  vs.  20  sec), 
the  error  appears  to  be  tolerable. 

One  may  have  noticed,  from  the  convergence  study,  that 
the  medium  model  is  nearly  as  accurate  as  the  fine  model 
for  the  non- ideal  case  of  full  tab  clamping  =  1  3/4") 
used  in  the  convergence  study.  Since  for  the  nonlinear 
analysis  the  boundary  conditions  are  more  ideal,  the  stress 
field  is  more  uniform  and  hence  the  fine  model  is  not 
required  to  give  a  good  overall  picture  of  the  specimen 


behavior.  Therefore,  the  medium  model  was  used  in  the 
nonlinear  analysis  of  the  specimen.  This  resulted  in 
approximately  a  40X  savings  in  cpu  time.  It  is  estimated 
that  a  nonlinear  run  using  the  fine  model  would  have  taken 
5000  cpu  sec  (computer  used  was  a  Cyber  170/845). 

The  above  results  will  be  discussed  in  depth  in 
chapter  5. 


CHAPTER  4 


& 


ii 


i 


■ 


trr 


EXPERIMENTATION 


4.1.  I ntroduct 1  on 


Two  sets  of  experiments  were  carried  out  In  the 
present  study.  The  first  set  of  experiments  were  simple 
tension  tests  to  obtain  the  basic  tensile  material 
properties  of  the  materials  used  to  fabricate  the  off-axis 
specimens.  The  second  set  of  experiments  were  off-axis 
tension  tests  using  a  new  fixture  which  could  apply  the 
load  with  the  ideal  boundary  conditions  determined  in 
chapter  3 . 


4.: 


Specimen  Fabrication 


The  material  system  used  in  this  study  was  AS4  3501-6 
graphite/epoxy.  The  material  was  supplied  by  Hercules 
Incorporated  in  the  form  of  1 2- in.  wide  prepreg  tape.  The 
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prepreg  material  was  used  to  fabricate  one  16-piy  panel, 

30  X  30  in.,  and  one  8-ply  panel,  12  x  12  in.  The  16-ply 
panel  had  a  lay-up  while  the  a-ply  panel  had  a 

(±45)22  layup.  Both  panels  were  cured  according  to 
manufacture' s  specifications.  After  curing,  both  panels 
were  subjected  to  C-scan  and  X-ray  inspection  for  flaws, 
before  being  cut  into  specimens.  The  Inspection  did  not 
reveal  significant  defects. 

A  single  panel  of  1/16  in.  thick  G-10  glass/epoxy 
(0/90  woven)  was  also  obtained  from  a  local  plastics  dealer 
for  tabbing  of  the  basic  material  properties  specimens  and 
the  off-axis  specimens. 

Specimens  for  basic  material  properties  were 
fabricated  from  the  above  panels  for  tension  tests  to 
determine  values  of  Ej,  E2*  ^12*  ''l2*  ’ '/e  specimens 

were  fabricated  for  each  fiber  direction.  Refer  to  Table 
4.1  for  the  tests  used  to  determine  the  various  material 
properties.  Note  that  the  commonly  used  ±45°  laminate  was 
used  for  determining  vs.  Yj2*  specimen  dimensions 

are  shown  in  Fig.  4.1  except  for  the  0°  graphite/epoxy 
specimens  which  had  a  reduced  width  of  1/2  in.  This 
reduction  was  required  because  the  1  in.  wide  0°  specimens 
could  not  be  loaded  to  failure  with  the  test  machine  used. 

Specimens  for  off-axis  tests  were  fabricated  from  the 
graphite/epoxy  panel.  Five  specimens  were  fabricated 
for  each  off-axis  angle.  The  specimen  dimensions  are  shown 
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4.2  Off-Axis  Tension  Specimens 


in  Fig.  4.2.  Note  that  the  0°  Fibers  in  the  tabs  are 
oriented  parallel  to  the  fibers  In  the  graph  1 te/epoxy.  As 
was  done  for  the  finite  element  model  in  section  3.2. 


4.3.  Basic  Material  Property  Tests 

All  specimens  were  instrumented  with  two  3-e1ement 
strain  gage  rosettes,  one  on  each  face,  at  the  center  of 
the  test  section  of  the  specimen.  An  Instron  Test  Machine, 
Floor  Model  TT-1115,  was  used  to  test  all  specimens  at 
ambient  environments.  The  maximum  load  for  this  machine  Is 
20,000  lbs.  The  crosshead  speed  was  0.05  in.  per  minute. 

For  each  test,  strain  was  read  from  all  six  strain 
gage  elements  and  load  was  read  from  the  Instron  Test 
Machine  load  cell  at  up  to  four  times  a  second  using  an 
analog  to  digital  converter  and  a  computer.  The  gage 
readings  from  the  two  faces  were  averaged  to  reduce  the 
effect  of  any  bending  or  torsion  that  might  have  been 
present  during  testing.  Strain  gage  element  numbering  and 
orientation  Is  shown  in  Fig.  4.3. 

To  obtain  a  desired  stress-strain  curve,  the  needed 
data  was  read,  manipulated,  and  plotted.  To  determine,  for 
Instance,  Oj  vs.  c^,  the  axial  stress  (Oj)  was  plotted 
versus  the  average  of  the  strains  from  the  #1  gages  (Ej) 
from  a  0°  tension  test.  All  five  experimental  Oj  vs.  Cj 
curves  were  then  placed  on  one  plot  and  an  average  curve 
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was  drawn  by  hand  through  the  ffve  experimental  curves. 


these  average  curves  are  shown  in  appendix  B  and  the 
average  engineering  elastic  constants  and  average  strengths 
are  tabulated  in  appendix  C.  The  engineering  elastic 
constants  (initial  values  of  Ej,  £2*  ^12*  ''12^  were 

used  in  the  linear  finite  element  analysis  and  the  average 
curves  were  represented  by  the  cubic  spline  functions  in 
the  nonlinear  finite  element  technique  described  earlier. 

4.4.  Off-Axis  Tension  Tests 


fitted  Into  the  side  plates.  Fig.  4.7  shows  how  the 
specimen  Is  placed  Into  the  clamp.  Note  the  knurled 
clamping  surface  of  the  clamp.  Fig.  4.8  shows  an  alignment 
fixture  used  to  maintain  alignment  of  the  specimen  while 
the  clamp  bolts  are  being  tightened.  The  fixture  Is  bolted 
down  so  that  the  spec  1 men  and  c 1  amp  rema In  f I xed  while  the 
clamp  bolts  are  being  torqued. 

Looking  at  Fig.  4.7,  one  can  see  that  the  Jaws  used  to 
clamp  the  specimen  tabs  are  of  fixed  dimensions.  This 


presents  a  slight  problem  when  one  desires  to  adjust  the 


amount  of  the  tab  that  Is  clamped  (D  )  and  still  allow  for 

c 

variations  In  the  pin  position  (Dp).  Ideally,  the  fixture 

would  allow  for  Independant  variations  of  0^  and  D^.  Since 

c  P 

with  the  present  fixture  this  Is  not  possible,  the  portion 


of  the  tab  that  was  not  to  be  clamped  was  milled  down 
approximately  0.015  in.  This  prevented  the  knurled  surface 
of  the  jaws  from  bitting  into  that  portion  of  the  tab. 


'Vi 


This  modification  to  the  specimen  is  not  recommended  as  a 
general  practice  In  off-axis  specimen  design.  Rather,  the 
Jaws  should  be  modified  Instead,  leaving  the  simple 


specimen  design  Intact. 

The  values  of  0  and  0  used  in  the  off-axis  test  were 
c  P 

1.0  in.  and  0.5  In.,  respectively.  The  choice  of  D^  =  1.0 


Is  obvious  when  one  looks  at  Fig.  3.9.  While  Dp  =  0.5  does 
not  correspond  to  the  absolute  minimum  value  of  for 
either  the  10°  or  14°  specimens  it  was  used,  nevertheless. 


for  two  reasons.  The  first  reason  Is  that  It  allowed  the 


technicians  to  apply  a  uniform  clamping  pressure  very 
conveniently  with  a  torque  wrench.  The  second  reason  Is 
that  It  allowed  for  all  four  clamping  bolts  to  carry  an 
equal  share  of  the  clamping  force.  With  the  specimen  moved 
forward  so  that  =  0.0,  nearly  all  of  the  clamping  force 

Is  carried  by  two  of  the  bolts.  Concern  over  the  specimen 

possibly  pulling  out  under  these  conditions  prompted  the 
change  to  0  =  0.5  in.  The  difference  in  that  resulted 

from  using  0^  =  0.5,  rather  than  the  absolute  optimum 

value,  was  1.04  7.  and  3.98  7.  for  the  10®  and  14®  specimens, 

respect  I ve 1 y . 

Various  experimental  stress-strain  curves  were 
obtained  from  the  off-axis  tests.  These  will  be  presented 
and  compared  with  the  stress-strain  plots  obtained  from  the 
nonlinear  analysis.  In  chapter  5. 

A  coating  of  photoelastic  plastic,  used  for  evaluating 
stress  uniformity,  was  applied  between  the  strain  gage  and 
the  tab  on  one  of  the  10®  and  one  of  the  14®  specimens. 

The  coating  was  cemented  to  the  specimen  surface  with  a 
uniformly  reflective  adhesive.  Note  that  the  coating  was 


only  applied  to  one  half  of  one  side  of  the  specimen.  The 
coating  specifications  were, 

t=  0.010  ±0.001  in. 


where 


t  =  thickness 

k  =  strain  optical  coefficient 
f  =  stress  optical  coefficient 
As  the  coated  specimens  were  loaded,  they  were 
illuminated  with  white  light  and  photographed  through  a 
polarizer  at  discrete  load  levels.  The  photograhs  will  be 
presented  and  discussed  in  chapter  5. 


CHAPTER  5 


RESULTS  AND  DISCUSSION 


5.1.  I ntroduct 1  on 


In  this  section,  results  from  the  analytical  and 
experimental  efforts  will  be  presented  and  discussed.  In 
general,  the  discussion  will  be  carried  out  in  terms  of  the 
goals  of  the  study.  Again,  these  are  1)  determination  of 
boundary  conditions  (D^  and  Dp)  that  produce  the  most 
uniform  state  of  stress  in  the  specimen  gage  length  and  2) 
limited  verification  of  the  aforementioned  failure 
cr i ter ia. 


5.2.  Results  of  the  Linear  Parametric  Study 


Referring  back  to  Figs.  3.7  through  3.9  one  can 
clearly  see  from  these  plots  that,  as  mentioned  earlier. 


attains  a  minimum  for  1  in.  of  clamping  and  the  pin 


67 


toward  the  rear  of  the  specimen  for  both  off-axis  angles. 

There  is.  however,  more  that  can  be  drawn  from  these  plots. 

First  is  the  fact  that,  regardless  of  the  amount  of 

clamping,  the  worst  possible  place  to  locate  the  pin  is  at 

the  inside  edge  (D  =  D  )  of  the  clamped  area.  In  all 

P  ^ 

cases  there  is  a  decrease  in  when  the  pin  is  moved  away 
from  the  edge  of  the  clamp.  For  the  the  10®  specimen,  when 
D  =  D_  =  1.750,  Ao  is  99.7  %  greater  than  the  minimum  Ao 

C  p  X  X 

for  that  amount  of  clamping.  One  can  also  see  that,  for 

all  values  of  0^  for  both  off-axis  angles,  once  the  pin  has 

been  moved  1/2  in.  away  from  the  clamp  edge  (D^  <  - 

1/2),  there  is  less  than  4  7.  subsequent  variation  in  Ao^. 

This  is  especially  interesting  in  light  of  the  suggestion 

by  several  researchers  (1,6,  &  10)  that  the  pin  should  be 

located  exactly  at  the  inside  edge  of  the  clamp. 

A  possible  reason  for  the  relatively  large  values  of 

Ao  when  0  =  i s  that  when  the  pin  is  at  this  location 

X  pc 

it  is  coincident  with  two  discontinuities  in  the  specimen, 
namely,  a  change  in  thickness  and  a  change  in  materials 
through  the  thickness.  These  discontinuities  have  the 
effect  of  restraining  nodal  displacements  in  the  local  area 
which,  in  turn,  causes  an  increase  in  stress  field  non¬ 
uniformity.  Since,  in  all  cases,  the  maximum  and  minimum 
axial  stresses  occured  in  elements  adjacent  to  the  tab, 
this  nonuniformity  has  a  direct  influence  on  Aa^.  Locating 
the  pin  at,  or  near,  these  discontinuities  places  yet 


another  restraint  on  the  local  displacements,  adding  to  the 
effect  of  the  discontinuities  and,  thus  further  perturbing 
the  uniformity  of  the  stress  field. 

One  might  ask  whether  this  phenomena  accurately 
represents  what  occurs  In  experimentation  since,  in  the 
test  fixture,  the  displacement  is  not  applied  at  a  discrete 
point.  Rather,  it  is  applied  to  the  entire  clamp  through 
the  pin  and  embedded  bearing.  Recall,  however,  that  in  the 
finite  element  model,  a  third  layer  of  steel  is  added  in 
the  clamped  region.  The  effect  of  this  relatively  rigid 
layer  is  to  effectively  apply  the  displacement  to  the 
entire  clamp,  just  as  is  done  in  the  experimental  fixture. 

The  reason  that  Ao^  decreases  when  the  amount  of 
clamping  is  reduced,  appears  to  be  mainly  due  to  an 
increase  in  specimen  flexibility  in  the  vicinity  of  the 
load  application.  The  portion  of  the  tab  that  remains 
undamped  serves  as  a  flexible  transition  between  the  rigid 
clamps  and  the  specimen  test  section.  The  distorted  finite 
element  models  shown  in  Fig.  5.1  and  5.2  show  very  clearly 
how  this  increase  in  flexibility  effects  Poisson  type 
contraction  near  the  clamps.  In  Fig.  5.1  one  can  see  that 
when  =  1.750,  Poisson  contraction  is  forced  to  take 
place  in  the  test  section,  giving  rise  to  an  increase  in 
stress  field  nonuniformity.  For  =  l.O  (Fig.  5.2),  on 
the  other  hand,  most  of  the  Poisson  contraction  and  hence, 
stress  field  non  uniformity,  takes  place  In  the  undamped 
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portion  of  the  tab,  leaving  the  test  section  less  affected 

by  the  rigid  clamp.  Stress  field  nonuniformities  in  the 

tab  are  relatively  unimportant  since,  due  to  the  increased 

strength  supplied  by  the  tab,  failure  is  not  likely  to 

occur  there.  Subsequent  experimentation  with  off-axis 

spec i mens  showed  that  fa i 1 ure  does  not  occur  i n  the  tabbed 

portion  of  this  specimen  design. 

Fig.  5.3  shows  graphically  how  stress  field  uniformity 

is  enhanced  by  this  increase  in  specimen  flexibility. 

Fig.  5.3a  shows  contours  of  when  =  Dp  =  1.750  while 

Fig.  5.3b  shows  the  contours  when  D  =1.0  and  D  =0.25 

c  p 

for  the  14°  specimen.  Visual  comparison  reveals  that  the 
stress  field  has  become  more  uniform  with  increased 
specimen  flexibility.  Notice  the  regions  of  maximum  and 
minimum  stress  indicated  in  these  figures. 

Finally,  looking  again  at  Figs.  3.7  and  3.8,  one  can 
see  that  sensitivity  of  to  variations  in  Dp  is  greater 
for  the  10°  specimen.  On  the  other  hand,  looking  at  Fig. 
3.9.  one  can  see  that  the  sensitivity  of  to  variations 
in  D^  is  greater  for  the  14°  specimen. 

5.3.  Qualitative  Assessment  of  Stress  Uniformity 

Qualitative  information  concerning  the  uniformity  of 
the  stress  field  in  the  gage  length  of  the  off-axis 
specimens  was  obtained  from  two  sources.  The  first  source. 


Region  of  Mininum  Stress 


xial  Stress,  c  ,  Contours.  Contour  Increment  =  100  psi 


was  photographs  of  specimens  coated  with  a  photoelastic 
plastic  material.  The  second  source  is  the  fractured 
spec i mens . 

From  the  theory  of  photoelastic  ana  lysis t  it  can  be 
shown  that  when  a  photoelastic  material  is  illuminated  with 
white-light,  the  shade  (or  combination  of  colors)  at  any 
point  in  the  material  will  be,  for  a  two  dimensional  model, 
a  unique  function  of  Oj  ~  ®2  difference  in  the 

principal  stresses).  If  a  coating  of  photoelastic  material 
is  applied  to  an  actual  structure,  made  uniformly 
reflective  by  use  of  a  reflective  paint  or  cement,  the 
colors  will  be  proportional  to  Oj  ~  °2  structure  as 

well  (Ref.  17). 

Since,  in  the  present  study,  the  state  of  stress  in 
the  gage  length  is  essentially  two  dimensional,  such  a 
method  of  ana  lysis  all ows  one  to  make  qua  1 i tat i ve 
judgements  as  to  the  uniformity  of  the  off-axis  specimen 
stress  field.  If  the  experimental  stress  field  is  in  fact 
uniform,  then  the  color  of  the  photoelastic  coating  should 
be  uniform. 

The  photographs  of  the  14°  specimens  are  shown  in  Fig. 
5.4  and  those  for  the  10°  specimens  are  shown  in  Fig.  5.5. 
The  load  levels  at  which  the  photographs  were  taken, 
accompanies  each  picture.  The  photoelastic  material  covers 
the  entire  width  of  the  specimen  between  one  of  the  tabs 
and  the  strain  gage.  In  each  picture,  the  tab  is  at 


the  top  and  the  strain  gage  is  at  the  bottom  (see  Fig. 

5.4a) . 

One  can  see  that  in  ail  of  the  photographs  the  stress 
field  is  essentially  uniform  between  the  tab  and  the  strain 
gage.  There  are  however  some  color  variations  that  require 
explanation.  The  first  is  the  oval  shaped  discontinuity  in 
the  color  field  in  the  lower  part  of  the  coated  area  in  all 
of  the  10°  specimen  pictures.  Upon  inspection  of  the 
specimen,  this  variation  was  found  to  be  due  to  poor 
bonding  between  the  plastic  and  the  reflective  adhesive. 

The  second  is  the  color  variation  present  near  the  bottom 
of  the  coated  area  in  all  of  the  14°  specimen  pictures. 

Upon  close  inspection  of  the  lower  boundary  of  the 
photoelastic  coating,  it  could  be  seen,  and  is  somewhat 
visible  in  the  pictures,  that  a  large  amount  of  adhesive 
had  squeezed  out  from  under  the  photoelastic  coating  such 
that  the  entire  free  edge  of  the  photoelastic  plastic  was 
embedded  in  the  adhesive.  This  appears  to  have  distorted 
the  stress  field  in  that  area  of  the  photoelastic  plastic. 
Finally,  one  can  see  that  there  are  bands  of  color 
appearing  parallel  to  the  fibers  for  both  specimens.  The 
exact  cause  for  this  variation  is  difficult  to  determine. 
However,  it  appears  to  be  a  manifestation  of  the  high 
stresses  occur ing  in  the  load  carrying  fibers.  Apparently 
the  bands  corresponding  to  fibers  that  are  closer  to  the 
surface  than  neighboring  fibers  and,  thus,  have  a  greater 
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Influence  on  the  color  field.  Notice  that  the  bands 
correspond  to  the  dominant  color  from  the  next  load  level 
(compare  Fig.  5.5e  with  5.5f).  This  indicates  that  the 
stresses  In  the  bands  are  higher  than  In  the  surrounding 
area.  Also  note  that  the  bands  occur  at  different  places 
for  the  two  specimens*  Indicating  that  the  variations  are 
likely  a  f unct I  on  of  mater  1 a 1  var 1 at  1 on  rather  than  a 
function  of  specimen  geometry. 

The  Important  observation  to  be  made  from  these 
pictures  Is  that  near  the  source  of  the  stress  field 
nonuniformity,  the  end  clamps  and  tabs,  there  Is  no 
apparent  region  of  stress  concentration  or  variation  for 
either  the  10®  or  14®  specimens.  As  one  looks  at  the  color 
field  Immediately  adjacent  to  the  tab,  one  sees  the  same 
colors  that  appear  1  In.  or  2  In.  away  from  tab. 

The  photographs  of  the  fractured  specimens  are  shown 
in  Figs.  5.6  and  5.7.  One  can  see  from  these  pictures  that 
for  all  but  one  of  the  10®  specimens,  failure  occured 
entirely  In  the  test  section  of  the  specimen.  Furthermore, 
for  eight  of  the  ten  specimens,  the  failure  surface  was  at 
least  3/8  In.  away  from  the  tabs.  This  I s  an  Indication 
that  the  stress  field  was  of  a  uniform  enough  nature  to 
prevent  failure  from  Initiating  from  a  particular  point  of 
stress  concentration  near  the  tabs. 

At  this  point  It  Is  Important  to  observe  that  the 
finite  element  technique  used  has  thus  far  proved  quite 
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accurate  In  predicting  the  behavior  of  the  off-axis  tension 
specimen.  It  has  been  used  to  determine  Ideal  values  of  D 

c 

and  Dp  which  have.  In  fact,  resulted  in  a  nearly  uniform 
test  section  stress  field.  This  uniformity  has  manifested 
Itself  quite  clearly  through  the  uniform  colors  of  the 
photoe 1 ast 1 c  coat  1 ng  and  cons  f  stent  test  sect  Ion  fall ures . 
With  this  In  mind  let  us  now  look  at  the  results  of  the 
non 1 1  near  ana  lysis  of  the  of f-ax 1 s  spec  1 men . 


5.4.  Results  of  the  Nonlinear  Analysis 


The  bulk  of  the  nonlinear  results  were  obtained  for 
comparison  with  experimental  data  and  will  be  discussed  as 
such  In  subsequent  sections.  It  Is  appropriate,  however-, 
to  observe  the  effect  of  material  nonlinearity  on  stress 
uniformity  at  this  F>olnt  since  It  cannot  be  directly 
measured  experimentally.  Fig.  5.8  shows  how  varied 
from  initial  1 oad I ng  to  fall ure .  In  this  plot  I  s 
plotted  as  a  percentage  of  midsection  axial  stress,  o^,  so 
that  the  effect  of  nonlinearity  can  be  more  readily 
observed.  One  can  see  that,  as  Indicated  In  Ref.  (11),  the 
uniformity  of  the  stress  field  Is  degraded  slightly  by  the 
material  nonlinearity.  For  the  10®  specimen  there  Is  a 
0.32  %  Increase  from  Initial  load  to  failure,  while  for  the 
14°  specimen  there  Is  a  0.76  %  Increase.  These  small 
increases  In  relative  to  the  midsection  axial  stress. 


.ilVm 


indicate  that  the  assumption  of  material  linearity  for  the 
design  of  the  off-axis  specimen  resulted  in  littlet  if  any, 
error  in  determining  D  and  D  . 

C 

5.5.  Comparison  of  Analytical  and  Experimental 
Stress-Strain  Responses 

Figs.  5.9  and  5.10  show  the  exper i menta 1  and 
analytical  stress-strain  curves  for  the  10®  and  14® 
specimens,  respect i v 1 ey.  The  experimental  curves  represent 
an  average  of  5  tests.  One  can  see  that  from  initial 
i oad i ng  up  to  approx i mate  1 y  60  %  of  the  ana  1 yt i ca 1  fail ure 
load,  correlation  is  quite  good  for  both  off-axis  angles. 
For  the  remainder  of  the  loading,  however,  the  correlation 
is  degraded.  Exact  errors  in  failure  stresses  and  strains 
v^i  1  1  be  compared  in  the  next  section  in  terms  of  the 
failure  criteria.  For  now,  a  reason  for  the  overall  lack 
of  correlation  Is  sought. 

Looking  at  Figs  5.9  and  5.10  one  can  see  that  the 
axial  stress  at  which  the  10®  responses  begin  to  diverge  is 
around  44,000  psi  while  for  the  14®  responses  divergence 
begins  at  around  32,000  psi.  Referring  back  to  eq.  2.4,  we 
can  determine  the  shear  stresses  in  the  fiber  direction, 
corresponding  to  these  axial  stresses. 


EScperimental-Analytical  Stress -Strain 
Curves  for  10“  Off-Axis 


For  the  10®  response, 

Tj2  =  44,000  X  cosdO®)  X  sin(10°)  =  7524.4  ps  i 
For  the  14®  response, 

^12  ~  32,000  X  cos (14°)  x  sln(14®)  =  7511.5  psi 

Thus  we  can  see  that  the  Initial  divergence  corresponds  to 
a  Fiber  direction  shear  stress  of,  nominally,  7500  psl. 

Now  referring  back  to  appendix  B,  to  the  experimental  shear 
stress  strain  curve  for  graphite/epoxy,  we  can  see  that  a 
shear  stress  of  7500  psi  corresponds,  approximately,  to  the 
point  where  the  one  dimensional  shear  stress-strain  curve 
becomes  highly  nonlinear.  Since  the  one-dimensional  normal 
and  transverse  stress-strain  curves  are  essentially  linear 
at  al 1  stress  levels,  and  since  the  correlation  in  the 
linear  range  is  quite  good,  it  follows  that  the  shear 
stress-strain  curve,  which  is  input  to  the  nonlinear  finite 
element  program,  is  cause  for  at  least  part  of  the  lack  of 
correlation  subsequent  to  the  two  axial  stresses  given 
above. 

Recalling  that  the  input  shear  stress-strain  curve  was 
obtained  experimentally  from  a  laminate  ((±45)25)  't;  seems 
that  it  would  be  safe  to  say  that,  for  this  case,  the  shear 
stress-strain  response  of  a  graphite/epoxy  laminate  is  not 
the  same  as  that  of  unidirectional  graphite/epoxy  in  the 
nonlinear  regime.  One  can  see  that  it  is  at  the  point 
where  the  (±45)25  laminate  shear  stress-strain  curve  begins 
to  exhibit  appreciable  nonlinearity,  that  the  responses 
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di verge. 

This  argument  Is  further  born  out  In  Fig  5.11  where 
experimental  shear  stress-strain  curves  are  shown  for  10° 
and  14°  off-axis  and  (±45 >25  laminate  along  with  analytical 
off-axis  curves.  One  can  readily  see  that  despite  the 
appjarent  uniformity  of  the  off-axis  specimen  stress  fieldt 
the  experimental  shear  response  of  the  off-axis  specimen 
does  not  continue  as  far  Into  the  nonlinear  regime  as  does 
the  (±45)22  laminate.  Furthermore,  beyond  approximately 
3000  psi  the  shear  modulus  of  the  off-axis  curve  Is  above 

that  of  the  (±45)22  l3*T»1nate.  In  short,  outside  of  the 
linear  regime,  the  shear  response  of  the  unidirectional 
off-axis  specimen  differs  somewhat  from  that  of  the  (±45)22 
laminate.  Notice  also  that  the  analytical  off-axis  shear 
stress-strain  curves  are  almost  exactly  the  same  as  the 
experimental  (±45)22  curve.  This  Is  an  indication  that  the 
finite  element  technique  was  operating  properly  given  the 
stress  strain-curves  Input  to  It. 

It  should  be  pointed  out  that  If,  In  the 
nonlinear  regime,  the  shear  stress-strain  response  of  the 
(±45) 2g  laminate  is,  in  fact,  different  from  that  of  the 
off-axis  specimen,  the  divergence  we  have  observed  is 
really  of  no  surprise.  In  choosing  to  analyze  the  10°  and 
14°  specimens,  we  have  chosen  off-axis  angles  where  the 
specimen  behaviors  are  essentially  controlled  by  their 
shear  characteristics.  Thus,  to  the  extent  that  the  shear 
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stress-strain  curve  input  to  the  finite  element  program 
is  inaccurate  in  the  nonlinear  regime,  the  predicted 
response  will  likewise  be  inaccurate. 

To  this  point,  we  have  observed  that  the  finite 
element  technique  itself  is  not  responsible  for  the  lack  of 
correlation  between  the  analytical  and  experimental 
responses.  We  have  seen,  however,  that  the  shear  stress- 
strain  response  of  the  (±45  ^2s  laminate  misrepresents  the 
off-axis  shear  response  and  consequently  causes  some  of  the 
lack  of  correlation.  The  extent  of  the  error  caused  by 
this  misrepresentation  and  the  error  caused  by  the  failure 
criteria  will  now  be  discussed. 

5.6.  Comparison  of  Analytical  and  Experimental  Failure 
Data 

One  can  see  from  Figs.  5,9  and  5.10  that  if  the 
failure  criteria  had  predicted  failure  at,  or  at  least 
closer  to,  the  failure  loads  of  the  experimental  specimens, 
the  errors  in  the  ultimate  strain  states  predicted  by 
analysis  would  have  been  fairly  small.  Table  5.1  shows  the 
errors  in  the  strain  states  corresponding  to  the 
experimental  failure  loads.  These  relatively  small  errors 
indicate  that  using  the  (±45)25  shear  stress-strain  curve 
is  not  the  primary  reason  for  the  overall  lack  of 
correlation.  Rather,  it  is  the  failure  criteria  that  is 


Let  us  look,  then,  to  see  how  an  earlier  (lower  stress 
level)  prediction  of  failure  might  have  been  affected  with 
the  given  criteria. 

Recalling  the  discussion  of  section  2.5,  the  m 
parameter  of  the  failure  criteria  was  taken  to  be  unity. 

It  was  not  clear,  however,  whether  this  was  in  fact  the 
best  choice  for  accurately  predicting  failure  of  the 
material  used  in  this  study.  To  observe  what  value  of  m 
might  have  been  a  better  choice,  let  us  plot  the  function, 

(Kj/Kj)'"+  (K2/K2)"’  +  <K6/K6)"’  =  1  (5.1) 

for  the  first  element  that  failed  (see  Fig.  5.12  for 
location  of  the  element),  for  m  *  1,  m  =  3/4,  and  m  =  1/2. 
These  plots  for  the  10®  and  14®  specimens  are  shown  in 
Figs.  5.13  and  5.14,  respectively.  Vertical  lines  are 
drawn  at  the  experimental  and  analytical  failure  stress 
levels.  One  can  see  that  for  m  =  1/2,  a  definite 
improvement  in  correlation  would  have  resulted.  For  the 
10°  specimen  the  error  in  failure  load  would  be  reduced  to 
0.0  7.  and  for  the  10®  specimen  the  error  would  be  reduced 
to  4.3  7,.  This  adjustment  of  m  i s  not  meant  to  be 
conclusive,  but  to  show  a  trend  In  the  effect  m  has  on  the 
predicted  failure  load. 

Let  us  now  look  at  a  second  parameter  of  the  failure 
criteria  that  may  influence  the  predicted  failure  load. 


Effect  of  Variation  of  m  Parameter 
on  Failure  Criteria  For  10**  Off-Axis  Specimen 
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Fig.  5-1^  Effect  of  Variation  of  m  Parameter  on 

Failure  Criteria  for  14°  Off-Axis  Specimens 
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Recalling  that  the  failure  of  the  10  and  14®  specimens  is 
primarily  in  shear,  the  first  two  ratios  of  eq.  5.1  are 
small  relative  to  the  third,  thus  the  value  of  Kg  (the  area 
under  the  one  dimensional  shear  stress-strain  curve),  has  a 
very  strong  influence  on  the  failure  criteria.  Referring 
to  Fig.  5.11,  one  can  see  that  the  areas  under  the  shear 
stress-strain  curves  for  the  experimental  off-axis 
specimens  (89.09  psi  and  74.82  psi  for  10®  and  14®, 
respectively)  are  much  less  than  the  area  under  the  shear 
stress-strain  curve  for  the  (t45)2g  laminate  (230.77  psi). 
Since  there  is  some  doubt  as  to  whether  the  (-^5)2g  shear 
stress-strain  curve  is  represent  I ve  of  the  unidirectional 
material,  let  us  again  plot  eq.  5.1  for,  m  *  1  and  Rg  equal 
to  the  area  under  the  experimental  off-axis  shear  stress- 
strain  curve.  These  plots  are  shown  in  Figs.  5.15  and 
5.16.  Note  that  the  actual  shear  stress-strain  curve  used 
by  the  finite  element  program  is  still  that  from  the 
(±45)25  laminate.  Only  the  value  of  has  been  changed. 
One  can  see  that  using  the  areas  under  the  off-axis  shear 
stress-strain  curves,  the  failure  load  would  be 
underestimated  in  both  cases.  Note  however  that  the 
predicted  failure  load  has  moved  in  the  direction  of  the 
experimental  failure  load.  Again,  using  the  off-axis 
values  of  Kg  is  not  a  conclusive  solution  to  the  problem, 
but  it  does  indicate  a  possible  method  of  improving  the 
resu I ts . 
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CHAPTER  6 


CONCLUSIONS 


It  Is  possible  to  separate  the  conclusions  drawn  from 

this  study  into  two  areas.  The  first  being  related  to  the 

determination  of  values  of  0  and  D_  which  produce  the  most 

c  p 

uniform  state  of  stress  In  an  off-axis  specimen. 

The  conclusions  drawn  from  this  first  area  are  the 
following: 

1)  The  length  to  width  ratios  necessary  for  accurate 
testing  of  off-axis  specimens  can  be  greatly  reduced. 

2)  Clamping  only  a  portion  of  the  tab  area  can  reduce 
the  differential  stress  in  the  gage  length. 

3)  Selectively  locating  the  pin  can  also  reduce  the 
differential  stress  in  the  gage  length. 

4 )  A  comb i nat I on  of  I dea 1  c 1  amp i ng  and  pin  1 ocat I on 
can  minimize  the  differential  stress  In  the  gage  length  to 
a  point  where  a  nearly  uniform  stress  can  be  produced  in  an 
economical  standard  geometry  off-axis  specimen. 


5)  Regardless  of  the  amount  of  clampingt  the  pin 
should  be  located  at  least  1/2  in.  behind  the  clamp  edge. 


6)  The  new  test  fixture  can  be  used  to  experimental ly 
produce  the  desired  uniform  state  of  stress  In  off-axis 
specimens  using  the  ideal  boundary  conditions. 

A  second  area  of  this  study  from  which  conclusions  can 
be  drawn  is  related  to  the  verification  of  the  strain 
energy  failure  criteria  and  nonlinear  analysis. 

The  conclusions  drawn  from  this  area  are  the 
f o 1 1 ow i ng : 

1)  Material  linearity  can  be  assumed  when  designing 
off-axis  composite  specimens  with  a  minimal  loss  in 
accuracy. 

2)  Using  the  (±45)25  laminate  shear  stress-strain 
curve  in  the  nonlinear  finite  element  analysis  resulted  in 
a  small  but  tolerable  error  in  correlation  between 
experimental  and  analytical  response  of  the  off-axis 

spec imen. 

3)  The  graphite/epoxy  off-axis  specimen  optimized  for 
shear  does  not  yield  as  much  of  the  shear  stress-strain 
curve  as  does  the  (i45)25  laminate. 

4)  Reducing  the  exponential  m  parameter  may  improve 
the  accuracy  of  the  strain  energy  failure  criteria. 

5)  Reducing  the  Kg  term  may  also  improve  the  accuracy 
of  the  strain  energy  failure  criteria. 
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Appendix  A 


Stiffness  and  Compliance  Matrices  (Ref.  (18)) 


Comp 1 1 ance  Coef f 1 c 1 ent s « 

S,i  =  1/Ej 

^12  "'’12/^1  ^  “''21/^2 

S22  =  I/E2 
^66  ^  ‘‘^^12 
^16  ®  ^26  ° 


Sjj  =  Sjjm  +  (2Sj2+Sgg)n  m  +  5^20 
Sj2  =  5^2  1  1  ■‘■^22~^66 

§22  =  j2+Sgg)n^m^  +  ^22^^ 

^16  "  ~  ^^^22'’^^12~^66^^  ^ 
§26  =  ^2Sj j -2Sj2-S6g)n  m  -  (2s22~2S j2-S65)nm 
^66  =  2(2Sjj+2S22-4Si2-Se6)^^"’^ 
where , 

n  =  s 1 n0  and  m  =  cos0 


Stiffness  Coefficients 


Qn  = 
Qi2  = 

Q22  = 
^66  = 
Qi6  = 

Qll  = 
Qi2  = 
O22  = 

Qi6  = 
^26  = 
^66  = 


El/(  1-''i2''21  ^ 

U2iEi/(  1-'>i2''21  ^ 

^2/  (  1  “w  j  2'*2 1  ^ 

Gi2 

Q26  =  0 

Qljm^  +  2(Qj2+2Qgg)n^m^  +  022^"^ 

(Ql  1 +022-4^66^ 

+  2(Qj2+2qgg)n^m^  +  Q22m'* 

^^1  ■*■  ^^12~*^22‘'’^^66^'^’ 
(Ql  1 -Ql2“2Qg6)fi^"’  ■*■  ^Qi2~Q22''’^*^66^'^'’ 
(Ql  l+q22-2Qi2-2Q66)"^"’^  + 


where 


Fig.  fl.2  Experimental  cr„  vs.  for  Graphite/Epoxy 


Experimental  cf  vs.  c.  for  Glass/Epoxy 


Fig.  B.6 


Experimental  T^p  vs.  "Yjp  for  Glass/Epoxy 
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Engineering  Elastic  Constants 
AS4  3501-6  Graphite/Epoxy 
E, 


I 


= 


12 


18.844 

X 

10^ 

ps  i 

1  .468 

X 

10^ 

psi 

0.910 

X 

10^ 

psi 

Vj2  =  0.280 


X,  Longitudinal  Tensile  Strength 

Y,  Transverse  Tensile  Strength 
S,  In-Plane  Shear  Strength 


285.6  ksi 
6. 194  ksi 
11.08  ksi 


G-10  Glass/Epoxy 


E,  = 


E.,  = 


12 


3.018 

X 

10^ 

psi 

2.890 

X 

10^ 

psi 

0.588 

X 

10^ 

ps  i 

'*12 

0. 

170 

X,  Longitudinal  Tensile  Strength 

Y,  Transverse  Tensile  Strength 
S,  In-Plane  Shear  Strength 


51.38  ksi 
35.75  ksi 
7.740  ksi 
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